Astronomy & Astrophysics manuscript no. 16168 


©ESO 2013 


January 1 1, 2013 





The solar magnetic field since 1700. 
II. Physical reconstruction of total, polar and open flux 

J. Jiang, R. H. Cameron, D. Schmitt, and M. Schiissler 



o 



Oh: 

I 

o 



> 
o 

O 



X 



Max-Planck-Institut fiir Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany 
Received / accepted 

ABSTRACT 

We have used semi-synthetic records of emerging sunspot groups based on sunspot number data as input for a surface flux transport 
model to reconstruct the evolution of the large-scale solar magnetic field and the open heliospheric flux from the year 1700 onward. 
The statistical properties of the semi-synthetic sunspot group records reflect those of the observed the Royal Greenwich Observatory 
photoheliographic results. These include correlations between the sunspot numbers and sunspot group latitudes, longitudes, areas and 
tilt angles. The reconstruction results for the total surface flux, the polar field, and the heliospheric open flux (determined by a current 
sheet source surface extrapolation) agree well with the available observational or empirically derived data and reconstructions. We 
confirm a significant positive correlation between the polar field during activity minimum periods and the strength of the subsequent 
sunspot cycle, which has implications for flux transport dynamo models for the solar cycle. Just prior to the Dalton minimum, at the 
end of the 18th century, a long cycle was followed by a weak cycle. We find that introducing a possibly 'lost' cycle between 1793 and 
1800 leads to a shift of the minimum of the open flux by 15 years which is inconsistent with the cosmogenic isotope record. 

Key words. Sun: heliosphere - Sun: surface magnetism - Sun: dynamo 



1. Introduction 

Information about the global properties of the solar magnetic 
field during past centuries provides important constraints for 
models of the Sun's global dynamo and is also relevant for stud- 
ies of the past terrestrial climate (Usoskin 2008). However, reg- 
ular synoptic direct observations of the large-scale solar sur- 
face field and space-based measurements of the interplane- 
tary magnetic field cover only a few decades. The aim of our 
work is to reconstruct the Sun's large-scale magnetic field and 
its heliospheric open flux from 1700 onwards on the basis of 
sunspot number records. The tool used for this reconstruction 
is a surface flux tr ansport model (SFTM, iBaumann et al 1 12004 
Ijiang et aDl2010bl) . In contrast to previous attempts at suc h re- 
constructions fe.g.. lWang et al.ll2005l: ISchriiver et al.ll2002h . the 
source input for our SFTM model consists of semi-synthetic 
records of sunspot groups whose properties (such as the dis- 
tributions of emergence latitudes, areas, and tilt angles) obey 
the empirical statis tical relationships w ith the activity cycle 
strength studied in Ijiang et aP (1201 Oal henceforth referred to 
as Paper I). The heliospheric open flux is extrapolated from 
the distribution of surfa ce flux using the current s heet source 
surface (CSSS) model (IZhao & H oeksema" 1 995 alibi) . The ob- 
tained open flux can be compared with recon structions based 
upon the geomagnetic gg index ([Russell 1975; Lockwood et alJ 
Il999t ISvalgaard & Cliveill2005l: IRouillard et ai.ll2007l) or upon 
cosmogenic radi onuclide data such as the concentration of 
^°Be in ice cores ('Bard et al."l997l'Caball ero-Lopez et al.ll2004 
iMcCrack en 2007: Steinhilber et al. 20 10). 

This paper is organized as follows. The SFTM, the treat- 
ment of its sources, and the CSSS extrap olation are d escribed in 
Section 2. Observations and results from ICameron et al. (201 0.) 
(hereafter referred to as CJSSIO) are used to validate the model 
in Section 3. The results including the polar field, the total and 



the open magnetic flux from 1700 onwards are presented in 
Section 4 and possible errors are assessed in Section 5. Our con- 
clusions are given in Section 6. 

2. Model description 

2.1. Solar surface flux transport model 

The evolution of the large scale magnetic field on the solar sur- 
face can be studied using two-dimensional flux transport mod- 
els (tPevore et al. 1984; Wang et al 1989; Mackay et al. 200^ 
ISchriiver etani2002l: iBaumann et al.ll2004l) . which describe the 
passive transport of the radial component of the magnetic field, 
B, under the eff ect of differential rotation, Q, meridio nal flow, v 
( Babcocklll961h . and turbulent surface diffusivity, tjh (iLeightonI 
11964 ). A slo w decay due to the fa ct that diffusion occurs in three 
dimensions (ISchriiver et al.l2002l) is modeled here in the manner 
described in Baumann et al.l (l2006h . introducing the parameter 

The governing equation of the surface flux transport model 
(SFTM) is 
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where S (A, (p, t) is the source term of the magnetic flux, which 
describes the emergence of bipolar magnetic regions as a func- 
tion of latitude A, longitude cp and time t. Dii], ) is a linear opera- 
tor describing the decay due to radial diffusion. The value of 77^ is 
discussed in Section|3] Concerning the horizontal diffusivity, rjH, 
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Fig. 1. 12-month running average of monthly Wolf (Rz, blue) 
and group {Re, red) sunspot numbers during 1700-2010. Rc in 
cycle 23 uses the values of Rz- 

we use the reference value in CJSSIO as t jh = 2 50 km^s ' . The 
profiles of both dif ferential rotation. Q ([Snodg rass 19831) and 
meridional flow, v dvan Ballegooiien et al.lll998l) are the same 
as in CJSSIO: Q = 13.38 - 2.30 sin^^ - 1 .62 sin"* /I (in degrees 
per day) and 

^^^^^|llsin(2.4^)ms-i for W < 75°, 
\0 otherwise. 

While we take the meridional flow velocity to be time- 
independent, latitudinal inflows towards the active regions are 
taken into account via a reduction of the tilt angle as described 
in the next subsection (see also CJSSIO; iJiang et aP 1201 Oct 
ICam eron & Schussler "2010^. For the initial field distribution, we 
follow van Ballegooijen etjl] ( Il998h and CJSSIO. The strength 
of the initial field is determined by the parameter Bo coiTespond- 
ing to the field strength at the poles. The values used are dis- 
cussed in SectionlH 

The solution of the linear Equation ([TJ may be expressed in 
terms of spherical harmonics 

oo / 

BiA, 0, = Z '^)- (3) 

l=\ m=-l 

The axial dipole (/ = \,m - 0) and equatorial dipole compo- 
nents (Z = 1, m = ±1) are considered in the presentation of our 
results. We use the same definitions as I Wang et al.l(l2005h to cal- 
culate the solar surface axial (Dax) and equatorial (Deq) dipole 
strengths. The polar field is defined as the average over a polar 
cap of 15° latitude extension. The total surface flux is obtained 
by integrating the unsigned magnetic field over the whole sur- 
face. 

2.2. Sources of magnetic flux 

Here we describe how we model the magnetic flux source term, 
S{A,(p,t), in Equation ([TJ. We use the semi-synthetic sunspot 
group records determined i n Paper I on the basis of either the 
group sunspot number, Rc dHovt & Schatten|[l998h or the Wolf 
sunspot number, Rz, the 12-month running averages of which 
are shown in Figure [T] from 1700 onwards. Each sunspot group 
in the semi-synthetic record is taken to represent a bipolar mag- 
netic region (BMR). The area of a BMR, being the sum of the 
umbral, penumbral and facular areas, is calculated from the area 
of the co rresponding sunspot gr oup using the empirical relation 
found bv lChapman et al.1 (Il997l) . The tilt angle of the BMR, a, 
is assumed to be 70% of that of the sunspot group. This factor 
(70%) was empirically determined in CJSSIO and reflects the 
effect of the localized inflows into active regions. 



In Paper I, the tendency of sunspots to appear in activ- 
ity nests was measured by considering which combination of 
uniformly distributed random emergence longitudes and com- 
pletely ordered longitudes match the degree of non-randomness 
in the sunspot record of the Royal Greenwich Observatory since 
1874. Here we follow the same approach and use both the semi- 
synthetic datasets with ordered and random emergence longi- 
tudes. How the resulting fields are combined is discussed in 
Section 3. 

The magnetic field distribution of each new BMR is given 
by B(A, (p) = B+{A, (/>) - B iA, (p) with 

7,±r, ^^ R /0.4A/?\' 2[l - cos /3^{A, cf>)] ^ 

B*(^,0) = BmaxI — ^1 exp{ — }, (4) 

where /3±{A, (p) are the heliocentric angles between the location 
of the sunspot group,(/l, 0), and the centers of each polarity, 
('^±, <P±)- ^ is the separation between the two polarities. We take 

1/2 

Af3 - 0.45A „ , where is the total area of the sunspot groups , 
and 6-4° (Ivan Ballegooiien et alj|1998t iBaumann et al 12001 . 
Bm-dx is the peak field strength of a BMR. Following CJSSIO, we 
take B,nax = 374 G. 

A/3 and the tilt angle, a, determine the latitudinal separation, 
Ay6 sin a, of the two polarities of a BMR. This quantity is impor- 
tant for the axial dipole moment of the global magnetic field and 
affects the polar field and the open flux during activity minima. 
The longitudinal separation, given by A^S cos q;(cos /t) ' , dom- 
inates the equatorial dipole moment and the open flux during 
activity maximum phases. 

Associated with the «* BMR in the semi-synthetic sunspot 
record is the time t„ when the BMR appears. The contribution 
to the source term S(A,<p,t) from the BMR is taken to be 
6{t - t„) {B*{A, 4>) - B {A, (f>)), where 6 is the Dirac delta func- 
tion. This corresponds to modeling the sunspot group emergence 
process as being instantaneous. 

2.3. Extrapolation into the lieliospliere 

The SFTM describes the evolution of the magnetic field on the 
Sun's surface. To obtain the heliospheric open flux we have to 
extrapolate the surface field outward. Since the often used po- 
tential field source surfa ce model does not well repr esent the 
Ulysses spacecraft data (ISchiissler & Baumannll2006l) . we use 
the b etter suited current sheet s ource surface (CSSS) extrapola- 
tion (IZhao & Hoeksemalll995allbl) . 

There are three parameters in the CSSS extrapolation, 
namely, one associated with the thickness of the current sheet, 
a, the radius of the cusp surface, R^s, and the radius of the source 
surface, R^^. The values for the three parameters are taken from 
the reference case in CJSSIO, i.e., a - Q.2Rq, R^s - 1.55Rq, 
Rss = lO.O/?0. The location of the cusp surface determines the 
amount of the open flux, which is calculated by the integration 
of unsigned magnetic field over the whole cusp surface. 

3. Comparison of the model for the time period 
from 1913 to 1986 

In CJSSIO, we used the SFTM with a source term based on 
the actually observed sunspot longitudes, latitudes, areas and 
cycle-averaged tilt angles to reconstruct the surface field and 
open flux for the period 1913-1986. The tim e evolution of the 
open flux derived from geomagnetic indices (iLockwoodlliool 
iLockwood et al.l l2009l) . and the reversal times of polar fields 
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Fig. 2. Open flux from CJSSIO (soHd black curve), from the 
semi-synthetic sunspot group record with random longitudes 
(dotted blue curve), and from a combination of semi-synthetic 
models with random and ordered longitudes Z?com = ( 1 ^ c)Z?ran + 
cBoi-d with c - 0.15 (red curve). The asterisks show the maxi- 
mum value for each cylce. Results using the semi-synthetic data 
are averages over 20 realizations. 

dMakarov et al.ll2003h were well reproduced. We therefore be- 
gin by comparing the results obtained using the current model 
based on sunspot numbers with those in CJSSIO. 

To do so we have to consider the tendency for sunspots to 
occur in activity nests at specific longitudes, because this affects 
equatorial dipole moment and the open flux around activity max- 
ima. We performed simulations with BMRs emerging at ran- 
domly distributed longitudes and simulations where all BMRs 
appear at longitudes 90° in the northern hemisphere and 270° in 
the southern hemisphere. Since both the SFTM and the CSSS 
extrapolation are linear, we can combine the fields from the two 
kinds of simulations using 

Bcom(/?cs, 0, Cf>, t) = (1 -c)B,,n(/?es, 6, cf>, f) + cBord(/?c,s, 0, cf>, t). (5) 

where c is a measure of the strength of the activity nesting. In 
Paper I we found that the value c = 0.15 is appropriate for a 
determination of the open flux. This is confirmed here. Figure 
|2] shows the evolution of the open flux from CJSSIO compared 
with cases with c = (purely random) and c - 0.15. The case 
with random longitudes generates too low open flux at solar 
maxima. The rms difference between the results from CJSSIO 
is 1.19 X 10'"* Wb when c = (the purely random case) com- 
pared to 1.09 X 10'^ Wb when c = 0.15. When we consider 
only the cycle maximum values of the open flux for each cy- 
cle, the difference is greater, with the rms deviation falling from 
1.93 X lO'"* Wb for c = to 1.00 x 10'^ Wb for c = 0.15. 

Next we compare our simulations during 1913 - 1986 (cy- 
cles 12-21) based on Rq and Rz with the results from CJSSIO 
and with the observed data. Note that all results shown are based 
on averages over 20 sets of random realizations of the semi- 
synthetic sunspot records. We first set the radial diffusivity to 
Tjr - 0. Figure [3ja) compares the time evolution of total flux 
from CJSSIO with the models based on Rq and Rz, respec- 
tively. The three curves almost overlap. For the period from the 
early 1970s onwards, the direct magnetic measurements from 
the Mount Wilson and Wilcox observatories are also shown. 
Figure [3lc) shows the evolution of the polar field for both hemi- 
spheres. Since the input of the BMRs emergences is randomly 
distributed on the two hemispheres, our reconstructed north and 
south polar fields are similar The reversal times (indicated by 



cyan v ertical lines) are similar to those given by iMakarov et al.l 
(l2003h . The largest differences between the models occur in cy- 
cle 19. Figure [3ld) shows the evolution of the modeled open 
flux in comparison with that inferred from the geomagnetic aa 
index (iLockwood et al.l [20091) . For the whole time period, the 
rms difference between the inferred values from the aa index 
and CJSSIO, Rz and Rq are 0.99, 0.92, 0.96 xlO''' Wb, respec- 
tively, correspon ding to about 14% of th e average value of 6.92 
xlO''^ Wb of the ILockwood et al.l ( l2009l) data. 

The time-latitude plot of the longitudinal averaged signed 
photospheric field (magnetic butterfly diagram) during this time 
period is shown in Figure |3lb) for the reconstruction based on 
Rz- As found in previous studies, the latitude separation of the 
two polarities leads to a net flux when azimuthally averaged. 
The advection of this flux to the poles reverses the polar fields 
each cycle. The regular polar field reversals and the anti-phase 
between the polar field and the low latitude field are well repro- 
duced. 

These results show that our SFTM with sources based on the 
sunspot numbers Rz or Rq is consistent with the reconstruction 
of CJSSIO, which used the recorded properties of the actually 
sunspot groups, and also compares well with the observed pho- 
tospheric and the heliospheric magnetic field. 

Be cause the sunspot number data are less reliable before 
1849 (IVaaueroll2007t ISvalgaard & Clivedl2010l) we have cho- 
sen to use a non-vanishing value for the radial diffusivity t;,.. 
This was found to be necessary because when rjr - Q is, used, 
the e-folding decay time due to rin alone is approximately 4000 
years. This timescale was obtained numerically, and is consider- 
ably longer than the simple estimate ttRq/j] because the merid- 
ional flow tends to keep the magnetic field at the two poles sep- 
arated. The long e-folding time means that the field at any given 
time is affected by errors in the sunspot numbers at all previous 
times. This is very undesirable because noisy data at the begin- 
ning of the dataset then contribute for the entire period covered 
by the simulations without significant damping. Introducing a 
weak radial diffusivity rj,- - 25 km^s~' (which lead s to a decay 
time about 20 yr, as found by iBaumann et al. I l2006h is aimed at 
keeping 77^ small while still being able to sensibly use the early 
data. Figure [3je) shows the polar field evolution with the weak 
radial diffusion included. Compared to Figure |3lc), the largest 
differences are on the order of 20%. The time evolution of the 
corresponding open fluxes is shown in Figure |3lf). During the 
minimum phases of some cycles, the model is now closer to the 
data inferred from the aa index (e.g. Rq model around 1955 and 
Rz model around 1976), while other deviate more strongly (e.g., 
Rz model around 1924 and 1936). For the whole time period, the 
rms deviations for the models based on Rz and Rq are 0.93 and 
0.96 in 10'"^ Wb, respectively. 

CJSS 10 found a strong correlation between the polar field of 
cycle n and the strength of the subsequent cycle « -1- 1 . Without 
radial diffusion, the correlation coefficients are 0.80 and 0.52 for 
the Rz and Re, cases respectively. In our case with a sample size 
of 7, a significance level of p - 0.05 corresponds to r = 0.74. For 
all cases, there is no correlation between the polar field around 
the activity minimum of cycle « and the strength of the same 
cycle. The introduction of radial diffusivity and the correspond- 
ing decay of the polar field decreases the correlation between 
the polar field and its subsequent cycle strength and somewhat 
increases the correlation between the polar field of cycle n and 
the strength of the same cycle. All these correlation coefficients 
do not exceed the level required for significance at the p - 0.05 
level. 
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Fig. 3. Comparison of the model results based on Rz (blue curves) and Rq (red curves) with the reconstruction of CJSSIO (black 
curves) for the period 1913-1986. Panels (a)-(d) show the case without radial diffusion, (a) Total surface flux. Square and plus 
symbols during the last 2 cycles denote the observations from the Wilcox and Mount Wilson solar observatories, respectively, (b) 
Time-latitude diagram of the longitudinally averaged signed surface magnetic field from the model based on Rz- (c) Polar field. 
Dashed and solid curves are for th e southern and norther n hemispheres, respectively. The cyan vertical lines indicate the inferred 
times of polar field reversals from iMakarpv et al.l (l2003b. (d) O pen heliospheric flux. The dash-dotted curve represents the result 
inferred from the geomagnetic aa index ( Lockwood et alj|2009h . (e)-(f): results with radial diffusivity 77, - 25 km-s (e) Polar 
field, (f) Open heliospheric flux. 



We conclude that the SFTM with input data based on sunspot 
numbers effectively describes the solar magnetic field since 
1913. Although the introduction of a non-vanishing radial diffu- 
sivity leads to some possibly undesired decay of the polar field, 
it decreases the error caused by the noisy early sunspot numbers. 
The photospheric and the heliospheric field are reasonably well 
reproduced. In the next section, we present the results back to 
1700. 



4. Results 

4.1. Time evolution of the reconstructed field 

Figures |4] and |5] show time series of various properties of the 
reconstructed solar magnetic field, based on Rq and Rz, respec- 
tively. Yearly values for the reconstructed open fluxes are given 
in Tables [T] and |2] (electronic version only). All results shown 
correspond to averages over 20 sets of random realizations of 
the semi-synthetic sunspot records. The evolution of the total 
unsigned flux (panels a), polar field (panels b), open flux (pan- 
els c), axial and equatorial dipole field strength (panels d) are 
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Fig. 4. Reconstruction based on Re during 1700-2010. (a) Total 
surface flux, (b) Polar field. Solid and dashed curves are for the 
northern and southern hemisphere, respectively, (c) Open helio- 
spheric flux. The red curve gives the reconstruction, the grey 
curve is the 11 -yr running ave rage. The open flux as inferred 
from the geomagnetic aa index (iLockwood et al.ll2009l) is repre- 
sented by the black curve, (d) Axial (solid curve) and equatorial 
(dashed curve) dipole field strength. 



shown. The first ~ 20 yrs of these series are still affected by the 
choice of the amplitude of the initial field, Bo- Since cycle -4 
(1700-1712) is very weak in the Rq data, we set Z?o = in this 
case. For Rz we used Bq = 3 G. 

The polar field displays regular reversals, except for the Rc 
case during the Dalton minimum, when a reversal appears to 
fail. The flux transported to the poles during the weak cycle 5 
just cancels the strong polar field generated by the strong and 
long cycle 4. On the other hand, a clear reversal occurs in the 
case based on Rz (cf. Figure |5]l. Owing to the uncertainties in 
the sunspot number record, it cannot be determined with any 
confidence whether the polar field actually reversed or not. The 
grey curves in Panels (c) of Figures|4]and|5]are the 1 1-yr running 
average of the open flux. The long-term trend is compared with 
other independent reconstructions in Section |43] 
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Fig. 5. Same as Figure|4] but for the reconstruction based on Rz- 



Panels (d) of Figures |4] and |5] show the time evolution of the 
axial dipole field, D^^, and the equatorial dipole field, Dgq. Z)eq is 
approximately in phase with the evolution of the sunspot number 
and the total flux. Since we define the polar field as the average 
over latitudes poleward of +75°, D^^ varies in phase with the 
polar field. Due to the persistent emergence of sunspot groups 
in different locations, the evolution of Deq is noisier than that of 
Dax. This is partly mitigated by the fact that we show the average 
over 20 semi-synthetic sunspot group record realizations. 



4.2. Correlations with sunspot numbers 

The relation between the maximum polar field at the end of a 
cycle and the cycle strength (highest sunspot number during the 
cycle) is shown in the left panel of Figure |6] The middle panel 
gives the relation between the cycle strength and the sum of the 
maximum polar field of the preceding and the same cycle. This 
represents the change in the polar field from one cycle to the 
next. The correlation is expected to be stronger in this case, be- 
cause the tilted BMRs first reverse the polar field of the old cy- 
cle and then built up that of the new cycle. The right panel of 
Figure |6] shows the relation between the maximum polar field at 
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Fig. 7. Relation between sunspot number and total surface flux 
during cycle maxima (upper panel) and during cycle minima 
(lower panel). Red symbols refer to the Re case and blue sym- 
bols to the Rz case. 



Fig. 8. Top: relation between maximum sunspot number and the 
maximum open flux during a cycle. Bottom; relation between 
minimum sunspot number and the minimum open flux during a 
cycle. 



the end of a cycle and the strength of the next cycle. All cor- 
relations are significant at the p < 0.05 level. From the results 
described in Section 3, we expect that for = the correlation 
between polar field and strength of the next cycle would be even 
stronger. Such a correlation is potentially relevant in connection 
with flux transport dynamo models (e.g., Chatterjee et al. 200-4 
iDikpati et al.ll2004l) . 

Figure [T] shows the relationship between the total unsigned 
flux and the sunspot number During activity maxima, many 
BMRs emerge and the polar fields are weak. Hence there is a 
strong correlation between the sunspot number and the total un- 
signed flux. During the minima fewer sunspots emerge, and the 
total surface flux has a stronger contribution from the polar re- 
gion. This leads to a weaker correlation between the total surface 
field and the sunspot numbers during solar minima. On longer 
timescales, the total flux is roughly proportional to the sunspot 
number 

Since the contribution to the field strength of the multipole 
of order I falls off with radius as r""^^', the lowest-order mul- 
tipoles dominate the a mplitude of the open heliospheric flux 
(IWang & Sheelevl2002h . The emergence of any given BMR may 
increase or decrease the open flux, depending on whether its 
dipole moment vector is orient ed s o as to reinforce or reduce 
that of the pre-existing field (Wang et al. 2000). The left panel 
of Figure [8] shows that the maximum value of open flux during 
a cycle and the maximum sunspot number are strongly corre- 



lated. This results from the large number of BMRs emerging 
during the solar maximum period together with activity nesting, 
which generate a strong equatorial dipole field. On the other 
hand, the correlation is much weaker between the minimum 
value of open flux over a cycle and the minimum sunspot num- 
ber This is because the open flux during minimum phases is 
dominated by the axial dipole field corresponding to the polar 
field that has accumulated over the cycle. Note that the maxima 
of t he open flux usually oc cur 1-2 yrs after sunspot maximum 
(cf. IWang & Sheelevll2002h while the minima of the open flux 
are almost in phase with the sunspot numbers. 

4.3. Comparison of open flux with other reconstructions 

Figure |9] shows the comparison of our 11 -yr running average 
of the reconstructed open flux since 1700 w i th res ults based on 
geomagnetic aa index by iLockwood et al.l (l2009l Loc09) and 
Svalgaard & Cliver (20K SCIO ), from the cosmogenic isotope 
'"Be data by McCracken[( l2007i McC07) and Steinhilber et al] 
(l2010l S ABM 10), and from othe r models bas ed on sunspot num- 
bers by IVieira & Solankil (1201 0[ VS 10) and IWang et al.l (l2005[ 
WLS05). Our reconstructed open flux refers to the unsigned ra- 
dial component of the magnetic field at the orbit of Earth. Since 
the values given in SCIO correspond to the heliospheric vector 
magnetic field amplitude, a factor 0.4 has been applied . This f ac- 
tor is close to the value used by ISvalgaard & Chvell (l2006l) to 
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Fig. 6. Left: relation between the maximum sunspot numbers of cycle « and the maximum of the polar field at the end of cycle n 
for the Rz case (blue symbols) and the Rq case (red symbols). Middle: relation between maximum sunspot numbers of cycle n and 
the sum of the maximum polar fields of cycles n-1 and n. Right: relation between the maximum polar field of cycle « and the 
maximum sunspot number of cycle « + 1. Correlation coefficients, r, and significance levels, p are indicated. 




Fig. 9. Comparis on of the re c onstru cted open flux with other reconstructions. Loc09: see iLockwoodl (l2003h and Lockwood et al. 
(|2009|) . WLS05: Wans et al.l (l2005l) : solid and dashed green are thei r models SI and S2, resp ectively. S CIO: ISvalgaard & Clivei 
hoiol) . A factor 0.4 used t o convert their data into open flux. VSIO: IVieira & Solankil (|2010|) . S ABM 10: ISteinhilber et al.l ( 1201 Ol 
PCA composite). McC07: iMcCrackenI (l2007h . All data are 11-yr running averages, except SABMIO which is a 25-yr running 
average. 



convert the magnitude of the vector field to the radial compo- 
nent. 



All reconstructions show a similar behavior in the second 
half of the time period covered. The values given by SABMIO 
are systematically higher than the others by approximately 10%. 
The minimum open flux over the period studied is about 1- 
2x10'"* Wb and the maximum value is about 8x10'"^ Wb. 



Before ~ 1800, there is a distinct difference between the re- 
sults based on '°Be data (McC07, SABMIO) and the other meth- 
ods. The values inferred from "^Be are particularly high from 
1720 to 1750, when the Rq data indicate rather low solar activ- 
ity. Although the Rz values are somewhat higher than Re during 
this time period, our reconstructed open flux is still about 30% 
weaker than that of McC07 and SABMIO. 



4.4. A 'lost cycle'? 

Since the years 1790-1794 are poorly covered by sunspot ob- 
servations, it has been suggested that the unusually long cycle 
number 4 (1784.7-1798.3) may actually consist of two shorter 
cycles, so that a weak 'lost' c ycle could be missin g in the ex- 
isting sunspot number records (Usos kin et al. 112001 1). There is an 
ongoing discussion on this topic tKrivova et al l200llAriil2008l 
lUsoskin et al.ll2009l) . so that it seems reasonable to investigate 
how adding a 'lost' cycle 4' would affect our reconstructions. 

Using Re data, we followed Usoskin et al. (2001) and 
changed the cycle minimum times to 1784.3 (cycle 4), 1793.1 
(cycle 4', the possible lost cycle) and 1799.8 (cycle 5). All other 
cycle minima times were left unchanged (for values, see Table 1 
of Paper I). We then generated new semi-synthetic sunspot group 
records according to the method of Paper I. In accordance with 
Hale's polarity laws, the magnetic polarity orientation of the 
sunspot groups in cycle 4' was taken to be the same as that of 
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Fig. 10. Time evolution of the reconstructed (southern) polar 
field during 1760-1880 after adding the 'lost' cycle indicated 
by the two vertical lines between 1793.1 and 1799.8 (red curve) 
and without the additional cycle (standard case, black curve). 
Both reconstructions are based on Tee- 



the original cycle 4 and the orientations for all cycles before cy- 
cle 4' were reversed. 

Figure [To] shows the comparison of the time evolution of the 
polar field between the cases with and without cycle 4'. The po- 
lar field generated during cycle 4 is only slightly diminished by 
the short and weak cycle 4' and no reversal takes place. The 
following cycle 5 has the same polarities as cycle 4 and thus 
strengthens the polar field even further Thus the addition of 
the extra cycle leads to a strong polar field being maintained 
throughout the whole period between about 1 790 and 1810. This 
is followed by a period of weak polar field between 1820 and 
1830. In contrast to these results, the standard case without extra 
cycle predicts weak polar field between 1805 and 1820. 

Comparing the reconstructed open heliospheric flux with the 
measured "'Be concentration in ice cores should, at least in qual- 
itative terms, give an indication which model is to be preferred. 
Figure [TT] shows the evolution of 11 -yr running average of the 
re constructed open fl ux in both cases together with the "'Be data 
of iBeer et al.l (Il990h . The introduction of the extra cycle shifts 
the minimum of the open flux from about 1810 to about 1825, 
which is clearly inconsistent with the "'Be data showing maxima 
(corresponding to low open flux) around 1810 and 1817. The 
same inconsistency arises also when comparing with the other 
reconstructions of the open flux shown in Figure |9] 

The strong unreversed polar field and the delay of the open 
flux minimum at the beginning of the 19th century by ~ 15 years, 
which result from the introduction of extra cycle 4', appear to be 
highly anomalous and seem to be ruled out by comparison with 
the "'Be record. Our reconstructions therefore argue against the 
existence of a 'lost' cycle 4'. 

5. Possible sources of error in the reconstruction 

Our reconstructions are based on the datasets of the sunspot 
numbers Rz and Re- Both of these have errors, partly due to 
sparse observations during the early period (Vaquero 20()7j and 
partly due to subjective definitions o f what to include by the dif - 
ferent observers over the centuries (ISvalgaard & Cliverll201 Oh . 
Errors in sunspot numbers obviously affect our reconstruction 
of the magnetic field. In the following, we analyze the effect of 
possible errors in the sunspot numbers on the reconstructed mag- 
netic field by considering a weak cycle (cycle 6, 1810.6-1823.3) 
and a strong cycle (cycle 19, 1954.3-1964.9) as extreme exam- 




1780 1790 1800 1810 1820 1830 1840 
time 



Fig. 11. Time evolution of the reconstructed open heliospheric 
flux (Il-yrs running average) during 1780-1840 with the 'lost' 
cycle (red curve) and without (black curve), in com parison to 
the measured '"Be concentration in ice cores (iBeer et al .■ ■199dl) . 
Note the inverted scale for the ""Be data. 



pies. We restrict this analysis to the reconstruction based on Rz, 
noting that the results for the Rc case are very similar. 

The effects of a change of the sunspot number data on the 
polar field and the open flux in cycle 6 (left panels) and cycle 
19 (right panels) are displayed in Figure [12] The upper pan- 
els illustrate a reduction of the sunspot number by 30% (dotted 
curves) and an increase by 30% (dashed curves), respectively, 
of the original values. The solid curves represent the unchanged 
data. The previous 5 yr before the cycle starts and the subsequent 
45 yr after the cycle ends are shown as well. 

To simplify the discussion, we denote cycle 6 and cycle 19 
as cycle «, cycles 7 and 20 as « + 1, and so on. The effect of the 
variation of sunspot number of cycle n on the following cycles 
has a periodicity of two cycles and an amplitude which decays 
with an e-folding time of ~ 20 years. To reveal the effects in 
more detail, we have broken the time period into 5 intervals. 
Part I corresponds to the rising phase of cycle n (from the start 
to the maximum). The polar field (shown in middle panels of 
Figure \T2\ is only weakly affected during this interval, during 
which the emerging BMRs reverse the polar field of the previ- 
ous cycle. Hence higher sunspot numbers cause the polar field to 
become weaker and to reverse earlier The variation of the open 
flux (lower panels of Figure [12} during this time period is domi- 
nated by the variation of equatorial dipole field which is directly 
connected to the sunspot number Due to the competing effects 
of the axial dipole field (polar field), the relative variation of the 
open flux is less than that of the sunspot number 

Period II corresponds to the decaying phases of cycle « (from 
the maximum to the end of cycle). Part III is the rising phase of 
cycle n + 1 . During these two periods, the polar field and the open 
flux are strongly affected. More sunspots produce a stronger po- 
lar field and vice versa. The relative variation of the open flux is 
similar to the relative change of the sunspot number during pe- 
riod II, but the effect is weaker for part III since during the rise to 
the next maximum the equatorial dipole moment of cycle « + 1, 
which is unchanged, begins to dominate. 

Period IV mainly includes the decaying phases of cycle n + l 
and rising phases of cycle n + 2. The polar field is still signif- 
icantly affected but less so than during the periods II and III, 
owing to the effect of radial diffusion. More sunspots in cycle « 
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2000 



2010 



Fig. 12. Effect of varying the sunspot number (upper panels) by an increase of 30% (dashed curves) and by a decrease of 30% 
(dotted curves) on the evolution of the polar field (middle panels, average absolute value of the northern and southern poles) and the 
open flux (lower panels). The weak cycle 6 (1810.6-1823.3, left panels) and the strong cycle 19 (1954.3-1964.9, right panels) of 
the Rz case are considered. The Roman numbers denote 5 time intervals which are affected by the variation of the sunspot numbers 
in diff'erent ways (see text). 



produce a stronger polar field during this time. In period V, the 
effect of a variation of the sunspot number in cycle n becomes 
less and less important owing to the decay caused by the radial 
diffiasion. 

In addition, the derivation of the sources for the SFTM de- 
pends on the spatial and temporal properties of BMR emer- 
gences. Because we derive the source in terms of correlations, 
the model has random components. Figure [T3] shows the mean 
values and standard deviations arising from this randomness, 
derived using 20 independent realizations of the semi-synthetic 
sunspot group records. Again the results for a weak cycle (cycle 
6) and a strong cycle (cycle 19) are shown. The standard devi- 
ations are small for the polar field, but more significant for the 



open flux during the maximum periods and reach up to 13% of 
the open flux. The average standard deviations of the polar field 
and the open flux over the whole simulation time series are about 
10% of the corresponding mean values. 

6. Conclusion 

We have provided a physical reconstruction of the large-scale so- 
lar magnetic field and the open heliospheric flux since 1700 with 
a surface flux transport model with sources based on sunspot 
number data and on the statistical properties of the sunspot 
groups in the RGO photoheliographic results. The model has 
been validated through comparison with reconstructions based 
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Fig. 13. The effect of the random component in generating the 
sources for the SFTM. The reconstructed polar field (upper pan- 
els) and open flux (lower panels) are shown for cycle 6 (left) and 
cycle 19 (right). The red vertical bars indicate the standard devi- 
ations corresponding to a set of twenty synthetic sunspot group 
records with different random numbers. 

on the actual sunspot group record and with directly measured 
or observationally inferred quantities. 

Our source term, S , is based on the semi-synthetic sunspot 
group record of Paper I. It in turn is based upon correlations be- 
tween the cycle amplitudes and cycle phase, and sunspot group 
areas, emergence latitudes and tilt angles. Hence while the sur- 
face evolution described by the SFTM model is linear, our model 
of the source term introduces nonlinearities into the reconstruc- 
tion. These nonlinearities partly explain why the reconstruction 
for the period from 1874 to 1976 has polar field reversals each 
cycle without the need to i nvoke variations in the meridional 
flow (IWang & Sheelevl2002h or extra terms which cause the field 
to decav dSchriiver etal This was also found in CJSSIO. 

The sunspot numbers are less reliable prior to 1874, and we 
needed to introduce a slow decay of the field to produce polar 
fields which reverse each cycle. 

The reconstructions considerably extend the basis for corre- 
lations studies, such as the relation between the polar field am- 
plitude during activity minima and the strengths of the preceding 
and subsequent cycles, with implications for dynamo models. 
Introducing a possibly 'lost' cycle at the end of the 18th cen- 
tury leads to a shift of the open flux minimum during the Dalton 
minimum which is incompatible with the '"Be record. 

Acknowledgements. We are grateful to Y.-M. Wang, L. E. A. Vieira, M. 
Lockwood, and F. Steinhilber for kindly providing their reconstructed open flux 
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Table 1. Reconstructed open flux based on Rc 



Year 



1700- 


-1709 


0.01 


0.08 


0.16 


0.72 


0.96 


0.95 


1.20 


1.05 


0.98 


0.94 


1710- 


-1719 


0.88 


0.76 


0.71 


0.69 


0.74 


0.97 


1.32 


2.29 


1.80 


3.42 


1720- 


-1729 


3.46 


4.00 


3.29 


3.11 


3.31 


2.88 


4.09 


4.56 


4.16 


3.86 


1730- 


-1739 


6.12 


4.88 


4.21 


3.57 


3.34 


3.55 


4.71 


3.55 


3.04 


4.82 


1740- 


-1749 


3.33 


4.51 


3.56 


3.30 


2.69 


2.53 


2.43 


3.21 


6.01 


5.37 


1750- 


-1759 


5.43 


5.81 


5.12 


5.04 


4.32 


3.99 


3.93 


4.26 


4.54 


5.40 


1760- 


-1769 


5.07 


5.47 


5.48 


5.46 


5.24 


4.83 


4.29 


4.30 


5.06 


5.65 


1770- 


-1779 


8.06 


8.07 


9.96 


6.59 


6.74 


5.43 


5.13 


5.65 


6.29 


6.57 


1780- 


-1789 


5.01 


5.39 


5.09 


3.92 


3.31 


3.20 


3.95 


4.77 


6.08 


5.99 


1790- 


-1799 


6.99 


7.18 


7.21 


7.68 


7.04 


6.57 


5.92 


5.50 


5.61 


5.20 


1800- 


-1809 


5.17 


6.35 


5.54 


3.68 


3.29 


3.01 


3.15 


1.46 


0.72 


0.28 


1810- 


-1819 


0.13 


0.13 


0.44 


1.27 


1.94 


2.59 


3.95 


4.15 


4.23 


4.62 


1820- 


-1829 


4.43 


4.00 


3.67 


3.42 


3.59 


3.73 


3.93 


3.93 


5.23 


4.84 


1830- 


-1839 


6.68 


7.10 


5.74 


5.07 


4.61 


5.23 


6.84 


7.31 


6.77 


5.66 


1840- 


-1849 


5.75 


5.11 


4.37 


4.17 


3.86 


3.86 


5.10 


4.49 


6.89 


7.32 


1850- 


-1859 


8.70 


7.32 


7.63 


7.66 


7.18 


6.69 


6.35 


6.32 


7.46 


7.08 


1860- 


-1869 


7.04 


6.59 


6.21 


5.67 


4.59 


4.07 


3.83 


3.83 


4.14 


5.55 


1870- 


-1879 


6.61 


5.96 


7.80 


8.09 


6.48 


6.08 


4.64 


4.52 


4.42 


4.11 


1880- 


-1889 


4.21 


5.06 


5.74 


5.44 


5.32 


5.46 


5.51 


4.17 


3.62 


3.49 


1890- 


-1899 


3.25 


4.21 


4.80 


5.53 


6.71 


7.15 


7.31 


7.13 


6.08 


5.43 


1900- 


-1909 


4.88 


4.71 


4.41 


4.50 


5.42 


5.66 


6.49 


6.21 


6.51 


6.17 


1910- 


-1919 


6.26 


5.35 


4.41 


4.02 


3.91 


5.20 


6.34 


7.07 


6.28 


8.24 


1920- 


-1929 


7.21 


6.23 


5.82 


5.22 


4.91 


5.35 


6.89 


6.20 


7.61 


6.27 


1930- 


-1939 


6.19 


5.43 


4.80 


4.58 


4.55 


4.67 


5.69 


8.38 


8.10 


9.69 


1940- 


-1949 


10.70 


8.99 


7.35 


6.68 


6.18 


6.17 


7.68 


9.42 


8.61 


7.85 


1950- 


-1959 


10.50 


7.63 


6.43 


5.97 


5.19 


5.24 


6.87 


8.39 


9.42 


10.31 


1960- 


-1969 


9.77 


8.59 


7.06 


5.92 


5.42 


5.14 


5.16 


6.37 


6.66 


7.29 


1970- 


-1979 


7.87 


7.48 


8.31 


9.03 


7.91 


6.71 


6.37 


6.28 


7.30 


8.82 


1980- 


-1989 


9.81 


8.14 


11.16 


9.30 


8.99 


6.79 


6.09 


6.17 


6.44 


8.30 


1990- 


-1999 


8.78 


8.92 


9.14 


8.12 


6.34 


5.83 


4.95 


4.71 


5.74 


6.17 


2000- 


-2009 


7.38 


7.89 


8.31 


8.30 


7.90 


6.73 


6.89 


6.12 


5.51 


5.22 



Table 2. Reconstructed open flux based on Rz 



Year 



1700- 


-1709 


1.79 


1.94 


2.59 


2.98 


4.30 


4.96 


5.42 


4.53 


4.24 


4.64 


1710- 


-1719 


3.87 


3.47 


3.30 


3.28 


3.33 


4.22 


4.52 


3.99 


5.34 


4.89 


1720- 


-1729 


5.19 


4.99 


5.18 


4.55 


4.45 


4.93 


6.86 


6.36 


8.01 


7.55 


1730- 


-1739 


6.85 


5.84 


5.14 


4.76 


4.58 


4.90 


5.80 


6.80 


6.32 


8.23 


1740- 


-1749 


7.70 


6.61 


5.77 


4.60 


4.66 


4.32 


4.56 


5.05 


5.81 


5.08 


1750- 


-1759 


6.30 


5.62 


6.34 


6.25 


5.17 


4.82 


5.00 


4.95 


5.01 


4.95 


1760- 


-1769 


4.80 


6.20 


6.74 


5.86 


7.19 


5.62 


5.37 


5.31 


6.47 


7.54 


1770- 


-1779 


7.34 


6.65 


6.02 


6.84 


5.24 


4.46 


4.40 


6.21 


8.05 


9.64 


1780- 


-1789 


8.58 


7.95 


7.23 


5.95 


5.24 


4.84 


5.96 


6.51 


6.96 


8.51 


1790- 


-1799 


9.82 


6.94 


7.30 


7.71 


6.97 


6.92 


6.50 


5.92 


5.66 


5.33 


1800- 


-1809 


5.58 


6.14 


5.69 


5.30 


4.95 


4.90 


4.54 


2.95 


2.48 


2.07 


1810- 


-1819 


1.81 


1.76 


1.93 


2.34 


2.11 


3.23 


4.74 


4.86 


4.66 


5.28 


1820- 


-1829 


4.91 


3.98 


3.69 


3.51 


3.63 


3.73 


4.49 


4.88 


4.67 


5.84 


1830- 


-1839 


7.83 


7.63 


6.16 


5.77 


5.43 


5.82 


7.17 


8.72 


8.26 


7.49 


1840- 


-1849 


8.57 


6.41 


5.51 


4.93 


4.48 


5.11 


5.17 


6.06 


8.37 


11.06 


1850- 


-1859 


9.33 


9.64 


8.72 


8.37 


7.98 


6.97 


6.36 


6.26 


7.20 


7.56 


1860- 


-1869 


7.22 


6.04 


5.55 


5.38 


5.67 


5.53 


4.69 


4.39 


4.63 


5.92 


1870- 


-1879 


6.62 


8.21 


9.37 


9.30 


7.73 


5.45 


4.74 


4.83 


4.23 


3.94 


1880- 


-1889 


4.44 


5.38 


5.56 


6.19 


5.42 


5.31 


5.24 


5.06 


4.18 


3.85 


1890- 


-1899 


3.71 


4.28 


5.64 


5.83 


5.46 


5.66 


7.11 


5.46 


5.85 


4.83 


1900- 


-1909 


4.62 


4.10 


4.00 


4.09 


4.84 


5.27 


4.78 


4.28 


4.39 


5.66 


1910- 


-1919 


5.55 


4.89 


4.36 


4.03 


3.92 


5.34 


4.60 


6.75 


6.79 


6.16 


1920- 


-1929 


6.53 


5.33 


4.96 


4.45 


4.38 


5.06 


5.83 


6.52 


5.70 


5.44 


1930- 


-1939 


6.53 


4.96 


4.59 


4.27 


4.15 


4.70 


5.37 


8.08 


7.62 


8.32 


1940- 


-1949 


8.78 


8.95 


8.31 


7.17 


6.65 


6.31 


7.22 


10.75 


9.25 


9.58 


1950- 


-1959 


9.36 


9.05 


8.03 


5.93 


5.30 


5.46 


7.30 


8.83 


10.47 


10.95 


1960- 


-1969 


10.78 


9.63 


7.63 


6.88 


6.10 


5.63 


5.55 


6.75 


7.94 


6.11 


1970- 


-1979 


8.04 


8.05 


8.31 


8.30 


7.47 


6.85 


6.38 


6.18 


7.22 


8.21 


1980- 


-1989 


10.27 


8.82 


10.72 


9.00 


7.92 


6.55 


5.95 


5.87 


6.44 


10.18 


1990- 


-1999 


9.07 


10.13 


9.40 


7.89 


6.71 


5.55 


5.00 


5.05 


5.79 


5.07 


2000- 


-2009 


7.54 


6.96 


10.36 


8.11 


7.19 


6.77 


6.60 


5.68 


5.37 


5.06 



